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Abstract 

'Macroscopic  elastic  moduli  governing  the  incremental  deformation  of  pressurized  lungs  are 
calculated  on  the  basis  of  a  micromechanical  study  of  an  individual  lung  element  (alveolus).  The 
spongy  lung  tissue  (parenchyma)  is  presumed  to  be  in  an  initial  state  of  isotropic  tension,  and  the 
lung  alveolus  is  idealized  as  a  pin-jointed  truss  in  the  shape  of  a  regular  dodecahedron.  The 
analysis  is  based  on  a  variational  statement  of  non-linear  structural  mechanics,  and  the  results  show 
how  the  moduli  depend  on  the  initial  tension  and  the  assumed  constitutive  behavior  of  the  idealized 
truss  members. 

f 

INTRODUCTION 

The  theory  of  lung  elasticity  has  received  much  study  during  the  last  two  decades,  with 
particular  attention  directed  to  the  establishment  of  overall  constitutive  relations  between  stress  and 
strain  in  lung  tissue.  The  most  common  approach  (e.g.  Lambert  and  Wilson,  1973;  Frankus  and 
Lee,  1974;  Fung  et  al.,  1978;  Vawter  et  ah,  1979;  Lanir,  1983;  Stamenovic  and  Wilson,  1985)  has 
been  to  contemplate  representative  substructural  lung  elements,  idealize  their  geometries,  specify 
their  elastic  properties  (often  via  strain-energy  functions),  assume  that  local  and  overall  strains  are 
identical,  and,  finally,  get  macroscopic  constitutive  relations  by  appropriate  averages  of  the  local 
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stresses.  The  earlier  work  of  Mead  et  al.  (1970)  and  Wilson  (1972)  was  more  in  line  with  the 
complementary  approach  of  loading  representative  elements  with  the  macroscopic  stress,  and 
computing  their  strains.  In  a  recent  paper,  Kimmel  et  al.  (1986)  returned  to  this  viewpoint,  and 
studied  the  elastic  behavior  of  pressurized  lungs  on  the  basis  of  a  model  in  which  a  typical 
microstructural  element  of  spongy  lung  tissue  (parenchyma)  was  idealized  as  a  regular,  pin-jointed, 
dodecahedral  space  truss  (Fig.  la).  An  initial  isotropic  tension  transmitted  to  the  parenchyma  by 
the  surrounding  pressurized  pleural  membrane  was  simulated  by  radial  forces  at  each  joint  that 
produced  identical  tensile  forces  in  the  truss  members.  The  model  was  then  used  in  a  calculation  of 
linear,  isotropic  relations  between  subsequent  small  increments  of  macroscopic  stress  and  strain  in 
the  parenchyma.  The  results  showed  the  dependence  of  the  elastic  constants  for  incremental 
deformation  on  the  initial  hydrostatic  tension,  and  on  the  tensile  constitutive  behavior  of  the 
idealized  truss  members.  The  elastic  moduli  found  were  in  fairly  good  agreement  with  those 
inferred  from  available  experimental  data,  but  the  calculations  produced  some  anomalous  results 
for  small  values  of  the  stiffness  of  the  element  members. 

In  this  paper  the  model  of  Kimmel  et  al.  is  reformulated,  and  a  rigorous  analysis,  based  on 
a  general  variational  principle  of  non-linear  structural  mechanics,  is  executed. 

MODEL  FORMULATION 

The  microscopic  air  sacs  (alveoli)  of  which  the  spongy  lung  parenchyma  is  composed  are 
roughly  polyhedral  in  shape,  and  bounded  by  thin  interfacial  membranes.  Besides  alveoli,  the 
parenchyma  is  also  pervaded  by  blood  vessels,  by  bronchioles  emanating  from  the  main  bronchial 
tubes,  by  fibers  surrounding  openings  in  the  alveolar  walls,  and  by  fibrous  tissues  extending  from 
the  surrounding  pleural  membrane  (Wilson  and  Bachofen,  1982).  We  will  pretend  that  the  elastic 
effects  of  these  extra  constituents  may  be  simulated  simply  by  an  appropriate  thickness  of  the  walls 
of  an  alveolar  conglomeration.  Since  each  alveolar  wall  is  shared  by  two  alveoli,  we  will  regard 
the  parenchyma  as  an  aggregate  of  individual  hollow  polyhedra,  the  faces  of  w'hich  have  half  the 
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actual  average  effective  wall  thickness.  Then,  in  accordance  with  standard  ideas  of  composite- 
material  theory,  we  define  the  macroscopic  strain  produced  by  a  prescribed  macroscopic  stress  as 
the  average  of  the  strains  experienced  by  the  individual  polyhedra  in  a  large  piece  of  parenchyma 
subjected  to  boundary  tractions  consistent  with  the  given  stress.  When  the  lung  is  inflated,  the 
alveolar  openings  allow  the  pressures  in  the  alveoli  to  equalize,  so  that  there  are  no  net  pressures 
across  the  individual  alveolar  walls,  but  a  transpulmonary  pressure  P0  is  applied  to  the  inside 

surface  of  the  pleura.  We  will  assume  that,  in  tum,  the  pleural  membrane  transmits  a  hydrostatic 
tension  P  <  P0  to  the  parenchyma,  that  this  initial  hydrostatic  stress  produces  a  purely  dilatational 

strain,  and  that  subsequent  stress  and  strain  increments  are  then  connected  by  isotropic  constitutive 
relations. 

In  this  paper,  we  will  bypass  the  conceptual  averaging  process  needed  to  determine  the 
macroscopic  strain  produced  by  a  given  stress,  and  calculate  instead  the  average  strain  experienced 
by  a  single  hollow  element  in  the  shape  of  a  regular  dodecahedron.  In  this  analysis,  each  face  of 
the  dodecahedron  will  be  subjected  to  a  traction  consistent  with  the  prescribed  uniform  stress  but, 
at  the  same  time,  the  surface  displacements  will  be  constrained  in  a  way  that  keeps  each  face  nearly 
flat.  To  be  precise  (see  Fig.  lb),  it  will  be  assumed  that  (i)  the  vector  displacement  at  the  center  of 
each  pentagonal  face  is  the  mean  of  the  displacements  at  its  five  comers;  and  that  (ii)  on  each  face 
the  displacement  varies  linearly  within  each  of  the  five  constituent  triangles  that  meet  at  the  center. 
Accordingly,  the  five  comer  displacements  suffice  to  specify  the  displacement  within  each 
pentagon.  It  is  easily  shown  that  for  small  displacements  satisfying  these  kinematic  assumptions 


-4- 


the  average  mean  curvature  of  each  face  vanishes.*  By  virtue  of  its  symmetry,  the  isolated 
dodecahedron,  when  loaded  and  constrained  in  this  way,  is  indeed  an  elastically  isotropic 
structure,  properly  reflecting  the  isotropy  of  the  random  aggregate  it  is  meant  to  represent. 

Following  Kimmel  et  al.,  we  make  a  further  simplification,  in  which  we  concentrate  all  of 
the  in-plane  load  carrying  capacity  of  the  dodecahedral  faces  into  elastic,  pin-jointed  rods  along  the 
edges,  and  adopt  for  these  members  the  constitutive  relation 


Here  T  and  L  arc  current  values  of  member  force  and  length,  dots  represent  time  derivatives,  and 
B  is  a  non-dimensional  stiffness  that  may  itself  depend  on  L . 

We  will  first  consider  a  uniform  hydrostatic  tension  P  applied  normally  to  the  faces  of  the 
dodecahedron,  and  calculate  the  incremental  bulk  modulus  K  defined  by 

K  =  P/(V/V)  (2) 

where  V  is  the  enclosed  volume.  Then,  with  the  hydrostatic  tension  kept  constant  and  normal  to 
the  deforming  faces,  we  will  apply  surface  traction  rates  corresponding  to  a  uniaxial  stress  a  (Fig. 
2)  and  calculate  the  Young’s  modulus 

E  =  c/e  (3) 

where  e  is  the  average  strain  rate  in  the  direction  of  o  suffered  by  the  dodecahedron.  The 
standard  connections  for  isotropic  bodies  will  then  give  Poisson's  ratio  v  and  the  shear  modulus 
G  . 


*  For  small  deflection  w(x,y)  normal  to  an  initially  flat  surface  of  area  A  and  boundary  C,  the  local 
mean  curvature  is 

-jLj  V2wdA  =-^f  |^ds 
A  C 

where  3w/9n  is  the  normal  derivative,  and  the  line  integral  around  the  perimeter  of  the  pentagon 
vanishes  when  the  conditions  (i),  (ii)  are  imposed  on  w. 

t 


BULK  MODULUS 


The  calculation  of  K  that  follows  is  essentially  the  same  as  that  of  Kimmel  et  al.  Suppose 
that  hydrostatic  tension  P  produces  identical  forces  T  in  the  30  members  of  the  dodecahedron. 
Then,  by  the  principle  of  virtual  work,  P  produces  rates-of-change  of  volume  and  length  that 
satisfy 


PV  =  30TL 

(4) 

Self-similar  growth  of  the  dodecahedron  due  to  P  implies  that 

V/V  =  3L/L 

(5) 

so  that 

PV  =  10TL 

(6) 

and 

P/P  +  V/V  =  T IT  +  L/L 

(7) 

Eliminating  L/L  and  T/T  from  (5),  (7),  and  (1)  gives 

P/P  =  v/v 

and  so  by  the  definition  (2) 

K/P  =  (B-2)/3 

(8) 

YOUNG  S  MODULUS 

With  the  hydrostatic  tension  P  maintained,  we  apply  vertical  tension  a  to  the  12  planar 
faces  of  the  dodecahedron  (Fig.  2  shows  plan  and  elevation  views  of  the  loaded  dodecahedron; 
geometrical  relations  are  listed  in  Appendix  A).  The  nominal  stress  O  is  specified  as  a  force  per 
unit  horizontal  projection  of  original  surface  area,  i.e.  the  area  corresponding  to  the  presence  of  P 
alone. 

While  a  is  applied  P  is  kept  normal  to  the  current  surface  orientations,  and  its 
magnitude  in  terms  of  force  per  unit  current  area  is  fixed.  To  begin  the  analysis  of  deformation 
due  to  o  ,  we  assert  the  principle  of  virtual  work 
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30 

I T  8A,  =  P8V+  Bn  (9) 

M 

where  Aj  represents  the  elongation  of  the  i*  member,  Tj  is  the  corresponding  force,  and  7t  is 
the  potential  of  the  applied  a  loading.  (The  procedure  that  follows  was  guided  by  an  earlier 

development  (Budiansky  1969)  of  equations  for  the  nonlinear  analysis  of  space  trusses.)  The 
elongation  Aj  is  defined  as  (Lj-L0) ,  where  Lj  is  the  current  member  length,  and  Lq  is  the 

reference  length  corresponding  to  the  initial  state  when  only  P  is  present.  The  virtual  changes  in 
volume  8V  and  the  virtual  elongations  8Aj  must  be  compatible  with  changes  in  joint 

displacements  at  the  comers  of  the  dodecahedron  and  with  the  previously  stated  kinematic 
assumptions  concerning  the  displacements  within  the  faces  of  the  dodecahedron. 

Since  there  are  only  three  distinct  member  responses  when  a  is  applied,  namely  (see  Fig. 

2)  those  of  (1)  AB  ,  (2)  BC  ,  and  (3)  CD  ,  we  can  replace  the  intemal-virtual-work  sum  in  Eq. 
(9)  by 

3 

loX  Tj  Aj 

i-l 

Also,  just  four  displacement  variables  suffice  to  specify  the  deformation.  These  are  (see  Figs.  2, 

3)  the  vertical  displacements  wj ,  w2  in  the  direction  of  the  unit  vector  k  of  the  joints  at  A  and  C  , 
respectively;  and  the  outward  displacements  vj  ,  v2  at  B,  C  ,  in  the  horizontal  direction 
nB  =  nc,  pointing  away  from  the  vertical  central  axis  of  the  dodecahedron.  The  joint  displacements 
elsewhere  follow  by  symmetry. 

The  T5A  terms  in  (9)  will  be  rewritten  in  terms  of  the  Lagrangian  strain  defined  by 

q  =  l(L2-L20)/L20 

=  A/Lq+  ■j  (A  /Lq  )  (10) 

and  the  modified  force  variable  Q  defined  as 

Q  e  [L0/(L0+  A)]T  s  (L0/L)T  (11) 

It  follows  that 

8q  =  8A/Lq+  A8A/l4)=  (L/l|)8A 


and 


T8A  =  LqQSti 


so  that  the  variational  statement  (9)  may  be  transformed  to 

10L0  ^  Qk  8t]  k  =  P&V  +  8rc  (12) 

Note  that  since  "H  =  lL/Lq2 ,  and  T/T  =  Q/Q  +  L/L ,  the  constitutive  relation  (1)  is  equivalent  to 

Q/Q  =  (B-1)L/L 

=  (B-1)(L0/L)2  rj  (13) 

The  advantage  of  r\  over  A  is  that  it  is  quadratic  instead  of  irrational  in  the  displacements. 

It  is  convenient  to  introduce  nondimensional  displacements  yj  (i  =  1,2, 3, 4)  defined  by 

typ  y2>y3  .yj  =  [wj.vj.wj.vjJ/Lq  (i4) 

Then,  for  the  three  member  types,  the  definition  (10)  leads  to 

Tl k  =  aki  yi  +  4  b kij  yi  yj  (k*  1,2,3)  (15) 

in  terms  of  coefficients  a^  and  b^j  =  bkjj  that  are  derived  and  listed  in  Appendix  A.  Here,  and 
hereafter,  we  use  summation  convention,  summing  from  1  to  4  over  repeated  subscripts. 

The  deformed  volume  V  of  the  dodecahedron  can  be  expressed  as  a  function  of  the  join! 
displacements  by  (Appendix  B) 

v  =  V1+  ci  yj  +  j  djjyj  yj  +  0(y3)]  (16) 

Here  V0  is  the  original  pressurized  volume,  and  djj  =  djj .  Finally,  in  Appendix  C  the  potential  of 

the  applied  tensile  stresses  is  derived  in  the  form 

it=oV0lfjyi]  (17) 

Substituting  Eqs.  ( 1 5)-(  1 7)  into  Eq.  (12)  and  dividing  by  PV0  ,  puts  the  variational 
statement  of  equilibrium  into  the  form 

fnP)  i  Qk(*ki+  bkijyj)5y.  -tci  +  djj  yj  +  0(y2)]5y,  =  (a/P)f,5y,  (18) 

VO/  k-l 


Here  we  used  Eq.  (6)  to  replace  PV(/(10Lq)  by  T0  . 


(19) 


For  o  — » 0 ,  yj  vanishes  and  -»  Tq  ,  so  that  (18)  becomes 

*ki  -Cil  5yi  =  0 


Since  (19)  holds  for  arbitrary  variations  6yj ,  we  should  have 

^aki-Cj  =  0  (i  =  1,2, 3.4) 


(20) 


k-1 


and  it  can  be  verified  that  these  four  equations  are  indeed  satisfied.  This  is  consistent  with  the  fact 
that  (19)  simply  asserts  the  principle  of  virtual  work  at  the  initial  state,  and  our  special  assumptions 
concerning  the  displacements  ensures  its  validity. 

Now  we  want  to  use  (18)  to  get  a  variational  statement  governing  the  displacement  rates 
yj  .  Differentiating  (18)  with  respect  to  time  (remembering  that  the  8yj's  are  variations  that  remain 

arbitrary)  gives 

3 

^  [Qk  (flki  +  bki  i  yj )+  Qkbkij  yj]  5yj 

-[djjyj  +  0(yy)]8yj  =  (0/P)fi6yj  (21) 


Now  let  o  -» 0,  so  that  yj  -» 0,  L  ->  ,  and  Qk  -*  T0 ;  then 

1-^-  X  ^  aki  +  X  bkij  yj  “  dij  yjl  Syj  =  (a/P) fj  8yi 


Lo  k-l 


k-l 


(22) 


Here,  by  Eq.  (13), 

Qk  =  T0(B-l)iik 

where,  by  differentiation  of  Eq.  (15),  and  letting  yj  -♦  0, 

■Hk  -  aki  y; 


(23) 


(24) 


Hence  (22)  implies  the  four  independent  equations 

[(B-l)qiJ+  Sjj  d,j  1  y  j  =  (o/P)f  j 


(i  =  1 ,2,3,4) 


(25) 


?K?3g 
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a 


s« = h  b‘ 


are  listed  in  Appendix  B.  The  dotted  quantities  in  Eqs.  (25)  now  refer  to  the  state  at  o  =  0  ,  and 


these  four  equations  constitute  a  symmetric  system  of  four  simultaneous,  linear  equations  for 
yj  (j  =  1,2, 3, 4) .  The  results  will,  of  course,  be  proportional  to  6/P .  Then,  to  define  the  average 


strain  rate  £ ,  we  write 


k  =  oVn  e 


and  by  Eq.  (17)  this  gives 


£=fiyi 


The  final  result  for  Young's  modulus,  given  by  Eq.  (3),  is 


E/P  = 


(fi  yj  > 


Thus  the  solution  of  Eqs.  (25)  will  provide  E/P  as  a  function  of  B  .  Then,  with  the  use  of  Eq. 


(8)  for  K/P  ,  the  standard  relations 


E  =  3K  (l-2v)  =  2G(l+v) 


among  Poisson's  ratio  v  ,  shear  modulus  G  ,  E  ,  and  K  for  isotropic  materials  may  be  used  to 


compute  v  and  G/P . 


RESULTS  AND  DISCUSSION 


The  results  for  E/P  and  G/P  are  plotted  versus  B  in  Fig.  4,  and  Poisson's  ratio  v  is 


shown  in  Fig.  5.  It  is  shown  in  Appendix  D  that  for  B  -> »  ,  the  limiting  values  become 


E/P  =  5(1+  V5  )/4 


=  4.045 


0/P  =  1.348 
v  =  1/2 


Ijljlj  j  .J  |  L  ||  .11.1  |  I 
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For  B  -» 2  ,  K  vanishes,  and  v  ->  -1  ,  because  isotropic  expansion  dominates  the  strain- 
rate  tensor  induced  by  0  ;  hence  E/P  ~  9K  ~  3(B-2)  for  B  2  .  On  the  other  hand,  the  shear 
modulus  does  not  vanish  at  B  =  2,  where  G/P  =  1.127.  But  for  B  <  2  ,  the  structure  is  unstable, 
and  the  values  for  G/P  would  not  be  observable. 

The  limiting  results  for  B  -♦ shown  in  (32)  are  somewhat  higher  than  those  of  Kimmel 
et  al.  (1986),  who  find  E/P  ->3.51.  But  for  B  less  than  ~5,  the  results  of  Kimmel  et  al.  deviate 
very  markedly  from  the  present  ones,  showing  anomalous  increases  in  both  E/P  and  G/P  as  B 
decreases.  This  is  probably  due  to  the  fact  that  the  various  engineering  approximations  made  by 
Kimmel  et  al.  do  not  cope  adequately  with  a  delicate  balance  of  elastic  responses  as  the  instability 
point  at  B  =  2  is  approached. 

To  corroborate  the  isotropy  of  our  model,  an  independent  computation  of  Poisson’s  ratio 
was  made  by  means  of  a  direct  calculation  of  the  average  strain-rate  tensor  in  the  dodecahedron 
(Appendix  E).  An  explicit  formula  (Eq.  (E8))  derived  for  v  in  terms  of  the  y/s  was  found  to 

provide  numerical  results  in  perfect  agreement  with  those  calculated  from  E/P  and  K/P  on  the  basis 
of  the  isotropy  relations  (31). 

We  also  studied  (Appendix  F)  the  effect  of  imposing  an  additional  constraint  on  the 
displacements,  keeping  each  face  of  the  dodecahedron  perfectly  flat  during  the  imposition  of  a  .  It 
turns  out  that  for  B  «>,  this  additional  constraint  makes  no  change  at  all  in  the  limiting  results 
(32).  The  reason  for  this  is  that  inextensionality  of  the  members  enforces  the  planarity  constraint 
automatically!  Furthermore,  the  imposition  of  planarity  produced  only  a  tiny  increase  -  less  than 
1/10%  -  in  E/P  over  the  full  range  2  <  P  <  ■»,  and  the  results  in  Figs.  4,  5  were  barely  changed. 
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Experiments  on  dog,  pig,  and  horse  lungs  by  Hajji  et  al.  (1979),  involving  local  external 
indentation  of  the  pleural  membrane  of  an  inflated  lung,  were  used  by  the  investigators  to  extract 
values  of  G/(l-v)  for  incremental  straining.  Based  on  the  somewhat  uncertain  value  of  P/P0  =  .7 
(Stamenovic  and  Wilson,  1985,  suggest  P/P0  -  3/4)  the  results  are  shown  in  Fig.  6,  where  the 
experimental  values  of  G/[P(l-v)]  found  are  compared  with  the  predictions  of  the  present  theory. 
Theory  and  expcr'^ent  are  generally  consistent  for  the  range  B  -  5-20,  but  it  should  be  noted  that 
near  the  lower  end  of  this  range  the  theoretical  values  of  v  (Fig.  5)  become  quite  low. 

CONCLUDING  REMARKS 

An  important  feature  of  the  present  theoretical  approach  is  that  by  focusing  on  incremental 
stress-strain  relations,  the  constitutive  description  of  the  lung  substructure  has  conveniently  been 
thrown  into  the  single  elastic  parameter  B,  which,  to  be  sure,  may  itself  depend  on  the  stress  level. 
Not  only  is  B  supposed  to  represent  the  effective  elasticity  of  a  fairly  complex  structure,  it  can  also 
be  presumed  to  incorporate  the  effects  of  surface  tension  in  the  alveolar  walls.  (That  is  to  say,  the 
basic  relation  given  by  Eq.  (1)  can  still  be  considered  to  apply,  operationally,  when  surface  tension 
is  present.) 

The  Finite  limits  approached  by  G/P  and  E/P  for  B  ->  °°  are  clearly  consequences  of  the  fact 
that  the  basic  pin-jointed  dodecahedron  is  a  statically  underdeterminate  mechanism  (rather  than  a 
structure)  in  the  absence  of  hydrostatic  tension.  By  the  same  token,  it  is  clear  that  the  use  of  a 
statically  determinate  (or  indeterminate)  truss  model  would  have  led  to  G/P,  E/P  -*  «>  for  B  -»  °°. 
Finally,  we  note  that  if,  in  contrast  to  the  present  approach,  overall  properties  were  calculated  via 
the  imposition  of  a  uniform  local  strain  field,  then  B  -> ««  would  imply  unbounded  stiffness  for  all 
models  —  underdeterminate  or  not 
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APPENDIX  A 
Lagrangian  Strains 


The  geometrical  parameters  of  the  regular  dodecahedron  shown  in  Fig.  2  satisfy  the  relations 

to  =  54* 
a  =  72- 

sin  to  =  cos  =  (l+V5)/4 
4  sin2  to  =  1  +  2  sin  to 
h  =  Pj  =  L0  (sec  to)/2 
P2=g  +  h  =  L0tanio 
r  =  L0  tan  to  sin  to 

and  the  dihedral  angle  Q  between  normals  to  adjacent  planes  is 

fl  =  arc  tan  2 


These  relations  are  useful  in  the  calculations  of  this  Appendix,  as  well  as  in  those  that  follow. 

Consider  a  bar  PQ  having  initial  length  l_o ,  and  let  tpq  represent  a  unit  vector  pointing  from 
P  to  Q.  Displacements  Up,  Uq  at  P  and  Q  produce  the  deformed  length  Lpg  satisfying 

LpQ  =  liLgt  +  U  Q~ U  p  II  (A1 ) 


The  Lagrangian  strain  (Eq.  (10))  is  then  given  by 


tIpq  = 


tpq-  (UQ-Up)  11UQ-Upll2 

Lo  +  2L2o 


(A2) 


We  will  now  apply  the  generic  relation  (A2)  to  find  the  strains  of  members  AB,  BC,  and 
CD  (Fig.  3).  The  pertinent  joint  displacements  are 


UA=  Wik  +  vjn  A 
UB  =  wj  k  +  vj  n  B 

Uc  =  w2  k  +  v2  n  c 
UD  =-w2k  +  v2  nD 


(A3) 


Since 


r|3  =  tIcd=  (w2/L0)csc  2co  +  -i(v2/L0)csc  2co+ 2(w2/L0)2+  (1-sin  co)(v2/L0)2  (A6) 

Rewritten  in  terms  of  the  y^s  (Eq.  (14)),  Eqs.  (A4)-(A6)  provide  the  following  results  for 
the  coefficients  ay  in  (15): 

0  4  cos  co  0  0 

sec  co  -esc  2co  -  sec  co  esc  2co 

0  0  2  esc  2co  esc  2co 

The  coefficients  b^j  vanish,  except  for 

b122=  4  cos2  co  (A8) 


The  matrix  of  coefficients  is 


n>2iii= 


Finally,  the  only  that  do  not  vanish  are 


*>333  =  4 


b344  =  2(1 -sin  to) 

The  coefficients  qjj  and  Sy  defined  by  Eqs.  (26)  and  (27)  are 


sec  2  to) 


=  v  4 


1 

-  (esc  co)/2 

-1 

(esc  to)/2 

-  (esc  co)/2 

3(csc20))/2 

(esc  to)/2 

-  (esc2  to)/4 

-1 

(esc  co)/2 

9  -  8  sin  to 

(csc4to)/l6 

(esc  u)/2 

-  (CSC2  £0)/4 

(esc4  to)/ 16 

(esc2  co)/2 

Isyl  - 


4  -  2  sin  co 

0 

-1 


0 

-1 

0 

3  -  2  sin  co 


(A  10) 


(All) 


APPENDIX  B 
Volume  Change 

Together  with  the  assumptions  made  concerning  the  piecewise  linear  displacement  patterns 
within  each  face,  the  comer  displacements  Wj,  vlf  w2,  v2  (Fig.  3)  determine  the  configuration  of 

the  dodecahedron  subjected  to  uniform  tension.  The  deformed  volume  can  be  written  as  (Fig.  lb) 

^  =  ^ABFGH  +  ^^AXB  +  ^BXC  +  ^CXD  1  (®  1 ) 

wherein  each  of  the  subscripted  terms  represents  the  volume  of  a  pyramid  with  its  apex  at  the 
center  of  the  dodecahedron,  and  the  subscript  denotes  the  base.  The  original  area  of  the  pentagon 


ABFGH  was 


A  0  =  (5/4)  Lq  tan  co 


and  the  altitude  of  the  corresponding  pyramid  was  (Fig.  2) 


r  =  L0  tan  co  sin  co 


Hence  the  deformed  volume  is 


VABFGH  =  5(L0tan  co  sin  co+  w,  )(L0+  2v,  cos  co)2  tan  co/12  (B3) 

To  facilitate  the  calculation  of  the  remaining  terms  in  (Bl),  we  use  the  orthogonal  unit  vectors 
(i  j,k)  shown  in  Fig.  2  to  write  the  original  position  vectors  of  A,  B,  C,  D,  E  and  X  as 
A ,  B  =  [±i  12+  j(tan  co)/2+k  tan  co  sin  co]L0 

C ,  E  =  [+i  sin  co+  j  tan  co  sin  co+  k  (esc  2  co)/4]L0 

(B4) 

D  =  [j  tan  co-k(csc  2co)/4]L0 


X  =  [2j+  k](tan  co (1  +  sin  co)]/5 


The  displacement  vectors  are 
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wA,wB  =±ivj  cosco+  jvj  sin  co+kwj 
wc,wE  =  Ti  v2  cos  to  +  jv2  sin  to+  k w2 
WD=  JV2 -kw2 

WX=  twA+  WB+  WC+  WD+  WE^ 

=  j[2vj  sin  co+  v2(l  +  2  sin  co)]/5  +  k  [2wj+w2]/5 


N^gX=  [A  +  wa]  x  [X+wx]  •  [B  +  wb]/6  (B6) 

and  similarly  for  the  remaining  two  pyramidal  volumes  in  (Bl).  Carrying  out  the  required 
straightforward  but  tedious  calculations  leads  to 

V=  VQ+  [Wj  (3  +  8  sin  co)/2  +  2vj  (1  +  sin  to)  +  w2  (-1  +  4  sin  co)/2 
+  v2  (1  +  6  sin  co)]Lq  tan  to 

+  [Vj  (2  +  sin  co)  +  v2  (2  +  3  sin  co)  +  2wj  (3  +  2  sin  to) 

+  Wj  v2(7  +  6  sin  to)  -  4w2v,  +  4v,  v2  sin  to 
+  w2  v2(5  +  2  sin  to)]  [  2LC  sin  to]/3 
+  [cubic  terms] 


where  the  original  volume  V0  is 


\{j  =  5Lq  tan2  co  sin  co 


The  coefficients  in  the  expansion  in  (16)  for  V/V0  in  terms  of  the  nondimensional  y's  are  then 


=  cot  co  (3  esc  +  8)/ 10 


=  2  cot  co  (esc  +  1  )/5 
=  cot  co  (-esc  +  4)/ 10 


=  cot  co  (esc  +  6)/5 
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APPENDIX  C 
Stress  Potential 


If  the  applied  a  stresses  in  the  k  direction  produce  the  surface  displacements  u,  then  the 


potential  of  the  stresses  is 


n  =  J  [n  •  okk  »u]dS 
S 


where  the  integration  is  over  the  undeformed  surface,  and  n  is  the  unit  normal  to  this  surface.  Our 
kinematic  assumptions  imply  that  (Fig.  lb)  the  face  ABFGH  remains  plane,  as  does  its  counterpan 
at  the  bottom  of  the  dodecahedron.  Hence  the  partial  potential  of  the  o-stresses  applied  to  these 


two  faces  is  simply 


PoA^] 


where  A0  is  the  original  area  of  the  pentagon. 

In  face  ABODE  (Fig.  lb),  the  midpoint  vector  displacement  wx  is  the  mean  of  the  fi\e 
comer  displacements,  the  average  over  each  triangle  of  the  linearly  varying  displacements  within 
each  triangle  is  the  mean  of  its  three  comer  displacements,  and  it  follows  that  the  average 
displacement  over  the  pentagon  is  also  wx.  Hence  the  potential  of  the  o-stresses  is  just 

(n-okk*wx)A0  (C3j 


But  (Fig.  3) 


and  the  contribution  (C3)  is 


wx  •  k  =  (2wj  +w2)/5 
n  •  k  =  cos  Q  =  1/V5 

oA0  cos  Q  (2wj  +w2 )/5 


Adding  10  such  contributions  to  (C2)  gives  the  total  potential 
n  =  2oA0  [(1  +  2  cos  Q)Wj  +  w2  cos  f2) 

=  2  oA0[(  1  +  2/V5)Wj  +  w2  /V5] 

Hence,  with  the  use  of  Eqs.  (B2)  and  (B7)  for  A0and  V0,  which  give 


%  >  .0 


v'lvi-.' •.Ivi-V-!'-'  AV.  ■>  ■ 


V  ».  \  *.  N. 
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APPENDIX  D 
Inextensional  Limit,  B  -> 


For  B  ->«>©,  the  members  become  inextensional,  and  therefore  the  strain  rates  given  by 
Eq.(24)  must  vanish: 

aklyi  =  0  (k  =  1,2,3)  (Dl) 

It  follows  that  for  B  -» <», 

yi  =  P,c  (D2) 

w'here,  in  accordance  with  the  a^  matrix  shown  in  (A7), 

Pj  =  ^5  ,  p2  =  0  ,  P3  =  1  ,  P4  =  -2  (D3) 


To  evaluate  the  constant  c,  substitute  (D2)  into  the  variational  statement  (22)  and  choose 
5y,=  Pr  This  gives 

c[s,j  ~  djj ] pj  Pj  =  (o/P)f .  p ■  (D4) 

and  then  the  limiting  value  of  E FP,  from  Eq.  (30),  becomes 

E/P=  (slJ-d1J)p,PJ/(fkPk)2  (D5) 

Substitution  of  the  values  of  s,r  dir  and  fk  given  in  Eqs.  (A  12),  (B9),  and  (C9),  respectively, 

provides  the  result  (after  considerable  reduction) 

E/P  =  5  sin  co 

=  5(1  +  V5)/4  (D6) 

and  this  was  corroborated  numerically  via  the  solution  of  Eqs.  (25)  for  very  large  B. 


APPENDIX  E 


Average  Strain  Rates 

The  average  strain-rate  tensor  E  in  a  body  V0  that  suffers  velocities  u  on  its  boundary  S  is 


e  =  J(un  +  nu  )dS  /(2V0) 


where  n  is  the  unit  normal  to  S.  If  we  let  £  =  k-E'k  denote  the  strain  rate  in  the  k  direction,  (El ) 
gives 

£=|j(k.u)(n.k)ds|/V0  (E2) 

which,  by  (Cl),  is  the  same  as  n/(V0o)  =  f^.  Hence,  direct  averaging  to  calculate  E/P  gives  the 
same  result  as  that  found  from  the  use  of  the  energy  definition  (28)  for  e.  However,  an 
independent  calculation  of  Poisson's  ratio  can  be  based  on  the  direct  calculation  of  the  mean 


transverse  strain  rate 


eT  =  ( i  •  e  •  i  +  j  •  e  •  j  )/2 


Since  n  •  i  =  0  on  ABCDE  (see  Fig.  2) 


et=1o|J  (n  •  j)(u  •  j)ds|  /(2V0) 

'  a  nnnr  * 


ABCDE 


n  ■  j  =  sin  £2  =  2/V5 


J  (u  •  j  )dS  =  (wx.  j)A0 


ABCDE 


=  [2vj  sin  0)  +  v2  (1  +  sin  d))]A0/5 


where  (B5)  has  been  used  for  wx.  Hence  (E4)  gives 


et  =  2 [ 2 y 2  sin  0)+  y4 C1  +  2  sin  co))[A0L0]/(V0  J5) 


=  [  y2(l  +  V5)  +  y4 (3  +  V5)][ cot  0)  esc  co]/(4V5) 
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APPENDIX  F 


Planarity  Constraint 

An  edge  view  of  the  pentagon  ABCDE  (Fig.  3)  is  shown  in  Fig.  7.  For  the  comers  of  the 

displaced  pentagon  A'B'CDEi'  to  remain  co-planar,  we  must  have 

h  +  w.  -  w,  g  +  2w, 

- 1 _ l _ = _ f _ l _  (Fl) 

h/2  +  (v2  -  Vj )  sin  0)  g/2  +  v2  (1  -  sin  to) 

In  terms  of  the  non-dimensional  displacements  y* ,  this  reduces  to  the  condition 

F  =  g.  yt+  i h  ij  y  i  yj  =  0  ^2) 

where 

gj  =  (sec  tu)/2 
g?  =  tan  to 

(F3) 

g3  =  -  (sec  co)/2  -  2  tan  to 
g4  =  -  (sec  to)/2 

hi4=h41=-l+2sin  co 

h  23  =  h  32  =  2(1  +  2  sin  co)  .  (F4) 

h  34  =  h  43  =  -  (1  +  6  sin  co) 

and  all  other  h^’s  vanish. 

To  impose  the  constraint  (F2),  we  introduce  the  scalar  Lagrangian  multiplier  A.  ,  and  add 

X5F=  X(g,+  hjj  y3 )5yi  (F5) 

to  the  left  side  of  the  variational  statement  (18).  It  follows  from  this  modified  variational  equation 
that  if  >^->  Xq  for  o->  0  ,  then 

X  ak,"ci  +  ^ogijsy.  =  0  (F6> 

k«  1  J 

and  therefore  (Eq.  (20))  =  0.  If  we  now  reproduce  the  stey  that  originally  led  to  Eq.  (25),  the 

corresponding  equations  for  the  constrained  problem  become 
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[(B-D^+^-dylyj  +  Xgi-Cd/PJfi  (i  =  1,2,3, 4)  (F7) 

where  X  refers  to  the  state  at  o  =  0.  These  four  equations,  together  with  the  differentiated 
constraint  relation 

g.  yi  =  0  (F8) 

that  holds  for  =  0,  constitute  a  symmetric  system  of  Five  linear  equations  for  y  1(  y2,  y3,  y4  and 
X,  and  their  numerical  solution  provided  the  results  discussed  in  the  text.  The  equations  of 
inextensionality  (Dl)  and  the  planarity  condition  (F8)  are  equivalent  to  the  homogeneous  relations 


" 

" 

0 

1 

0 

0 

’  l 

2  sin  co 

-1 

-  2  sin  co 

1 

*2 

0 

0 

2 

1 

Y3 

1 

2  sin  co 

-  1  -  4  sin  co 

-1 

y4 

_ 

_ 

_ 

(F9) 


Remarkably,  the  matrix  of  coefficients  is  singular,  so  that  the  limiting  solution  (D2),  (D3) 
for  y,  in  the  limit  B->  °°  remains  valid,  as  do  the  limiting  results  for  the  moduli.  Thus,  planarity,  at 

least  for  small  displacements,  is  a  consequence  of  inextensionality. 


